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We  have  carried  out  high  quality  three-dimensional  ab-initio  quantum  chemi¬ 
cal  calculations  for  these  oximes  as  a  function  of  geometry  (CNOH  dihedral 
angles).  This  gives  us  calculated  conformational  energy  profiles.  All  pos¬ 
sible  conformations  within  a  certain  number  of  kcal  (on  the  order  of  12 
kcal/mole)  from  the  absolute  energy  minimum  may  be  possible  to  achieve  in 
solution  or  under  the  influence  of  an  active  site,  and  thus  may  be  significant 
pharmacol ogi cal ly .  The  equilibrium  geometry  of  a  species  is  where  the  total 
energy  is  at  a  minimum.  Our  results  of  quantum  chemical  calculated  total  energy 
vs.  dihedral  angle  for  syn  and  anti  2-PA  showed  that  there  were  multiple  maxima 
and  minima  in  these  compounds--thus ,  many  points  must  be  computed.  [For  these 
ab-initio  quantum  chemical  calculations,  we  are  using  our  own  fast  ab-initio 
programs  which  also  incorporate  optional  desirable  computational  strategies: 
MODPOT  -  effective  core  model  potentials  which  allow  calculation  of  the  valence 
electrons  only  explicitly,  yet  accurately;  and  VRDDO  -  (variable  retention  of 
diatomic  differential  overlap),  an  efficient  charge  conserving  integral 
prescreening  evaluation  which  tells  reliably  whether  an  integral  will  be  laroe 
enough  to  be  significant  or  so  small  that  it  should  not  be  calculated--which 
cut  down  the  computer  time  by  at  least  an  order  of  magnitude  even  for  a  com¬ 
pound  with  one  aromatic  ring  and  a  side  group  like  2-PA  and  considerably  more 
for  larger  molecules;  and  MERGE,  to  save  invariant  skeletal  integrals  and  reuse 
them.  This  MERGE  technique  is  especially  effective  when  the  geometry  changes 
and  also  saves  orders  of  magnitude  of  computer  time  when  series  of  related 
compounds  are  being  studied. 

As  we  calculated  the  quantum  chemical  wave  functions  of  those  oximes,  we 
generated  electrostatic  molecular  potential  contour  (EMPC)  maps  of  each 
molecule. 

So  far,  we  have  carried  out  such  calculations  for  syn  and  anti  2-PA  (2- 
pyridine  aldoxime),  syn  and  anti  4-PA  (4-pyridine  aldoxime),  syn  and  anti  2-PAM 
(2-pyridinium  aldoxime),  syn  and  anti  4-PAM  (4-pyridinium  aldoxime),  HI-6  (a 
bis-pyridium  monooxime)  and  imidazole,  protonated  imidazole,  and  BHH-1  (a 
methyl  imidazolium  aldoxime). 

We  have  also  carried  out  similar  ab-initio  quantum  chemical  calculations 
for  and  generated  EMPC  maps  around  deprotonated  2-PAM  (the  oximate  anion  sug¬ 
gested  as  possibly  being  the  active  species). 

The  three  dimensional  isopotential  contour  EMPC  map  for  the  deprotonated  2- 
PAM  showed  vividly  a  most  significant  finding.  The  three  dimensional  isopoten¬ 
tial  contour  EMPC  map  of  deprotonated  2-PAM  indicated  clearly  that  for  this 
species  there  is  a  strongly  positive  region  and-a  strongly  negative  region. 

This  has  .significant  implications  for  several  aspects  of  the  use  of  oxime 
reactivators. 

In  addition  to  our  theoretical  work,  we  have,  been  carrying  out  experimental 
determinations  of  the  1 ipophi 1 icities  (their  tendency  to  partition  from  aqueous 
to  lipid  phases)  of  these  oxime  reactivators.  We  have  carried  out  these 
1 ipophi 1 icity  determinations,  first  with  our  microelectrometr ic  technique.  It 
turns  out  that  neither  these  quaternary  oximes  nor  their  deprotonated  species 
are  very  lipophilic.  However,  since  there  have  been  reports  that  some  oxime 
reactivators  are  effective  at  low  concentrations,  we  have  set  up  and  are  also 
using  a  more  sensitive  UV  spectrophotometr ic  procedure  to  establish  quantita¬ 
tively  and  directly  just  how  much  of  these  compounds  is  partitioning  into  the 
lipid  phase. 

We  have  studied  the  following  compounds:  Microelectrometric  technique:  2- 
PAM,  HI-6,  TMB-4,  T0X0G0NIN,  BHH-1,  BHH-3,  BHH-4 ,  BHH-5 ,  BHH-8 ;  UV 
spectrophotometr ic  technique:  2-PAM,  BHH-5,  BHH-13,  BHH-14,  BHH-15. 


SUMMARY 


We  have  carried  out  high  quality  three-dimensional  ab-initio  quantum 
chemical  calculations  for  these  oximes  as  a  function  of  geometry  (CNOH 
dihedral  angles).  This  gives  us  calculated  conformational  energy  profiles. 
All  possible  conformations  within  a  certain  number  of  kcal  (on  the  order  of 
12  kcal/mole)  from  the  absolute  energy  minimum  may  be  possible  to  achieve  in 
solution  or  under  the  influence  of  an  active  site,  and  this  may  be  sig¬ 
nificant  pharmacologically.  The  equilibrium  geometry  of  a  species  is  where 
the  total  energy  is  at  a  minimum.  Our  results  of  quantum  chemical  calcu¬ 
lated  total  energy  vs.  dihedral  angle  for  syn  and  anti  2-PA  showed  that 
there  were  multiple  maxima  and  minima  in  these  compounds--thus,  many  points 
must  be  computed.  [For  these  ab-initio  quantum  chemical  calculations,  we 
are  using  our  own  fast  ab-initio  programs  which  also  incorporate  optional 
desirable  computational  strategies:  M00P0T  -  effective  core  model  poten¬ 
tials  which  allow  calculation  of  the  valence  electrons  only  explicitly,  yet 
accurately;  and  VRDDO  -  (variable  retention  of  diatomic  differential 
overlap),  an  efficient  charge  conserving  integral  prescreening  evaluation 
which  tells  reliably  whether  an  integral  will  be  large  enough  to  be  sig¬ 
nificant  or  so  small  that  it  should  not  be  calculated--which  cut  down  the 
computer  time  by  at  least  an  order  of  magnitude  for  a  compound  with  one 
aromatic  ring  and  a  side  group  like  2-PA;  and  MERGE,  to  save  invariant 
skeletal  integrals  and  reuse  them.  This  MERGE  technique  is  especially 
effective  when  the  geometry  changes  and  also  saves  orders  of  magnitude  of 
computer  time  when  series  of  related  compounds  are  being  studied. 

As  we  calculated  the  quantum  chemical  wave  functions  of  those  oximes, 
we  generated  electrostatic  molecular  potential  contour  (EMPC)  maps  of  each 
molecule. 
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So  far,  we  have  carried  out  such  calculations  for  syn  and  anti  2-PA  (2- 
pyridine  aldoxime),  syn  and  anti  4-PA  (4-pyridine  aldoxime),  syn  and  anti  2- 
PAM  (2-pyridinium  aldoxime),  syn  and  anti  4-PAM  (4-pyridinium  aldoxime),  HI- 
6  (a  bis-pyridium  monoxime)  and  imidazole,  protonated  imidazole,  and  BHH-1 
(a  methyl  imidazolium  aldoxime). 

We  have  also  carried  out  similar  ab-initio  quantum  chemical  calcula¬ 
tions  for  and  generated  EMPC  maps  around  deprotonated  2--PAM  (the  oximate 
anion  suggested  as  possibly  being  the  active  species). 

The  three  dimensional  isopotential  contour  EMPC  map  for  the 
deprotonated  2-PAM  showed  vividly  a  most  significant  finding.  The  three 
dimensional  isopotential  contour  EMPC  map  of  deprotonated  2-PAM  indicated 
clearly  for  this  species  that  there  is  a  strongly  positive  region  and  a 
strongly  negative  region.  This  has  significant  implications  for  several 
aspects  of  the  use  of  oxime  reactivators. 

In  addition  to  our  theoretical  work,  we  have  been  carrying  out  ex¬ 
perimental  determinations  of  the  1 ipophil icities  of  these  oxime 
reactivators.  We  have  carried  out  these  1 ipophi 1 icity  determinations,  first 
with  our  microelectrometr ic  technique.  It  turns  out  that  neither  these 
quaternary  oximes  nor  their  deprotonated  species  are  very  lipophilic. 
However,  since  there  have  been  reports  that  some  oxime  reactivators  are 


effective  at  low  concentrations,  we  have  set  up  and  are  also  using  a  more 
sensitive  UV  spectrophotometr ic  procedure  to  establish  quantitatively  and 
directly  just  how  much  of  these  compounds  is  partitioning  into  the  lipid 
phase. 

We  have  studied  the  following  compounds:  Microelectrometric  technique: 
2 -PAM,  HI-6,  TMB-4,  T0X0G0NIN,  BHH-1 ,  BHH-3,  BHH-4,  BHH-5,  BHH-8,  BHH-18;  UV 
spectrophotometr ic  technique:  2-PAM,  BHH-5,  BHH-13,  BHH-14,  BHH-15. 


FOREWORD 


Citations  of  commercial  organizations  and  trade  names 
in  this  report  do  not  constitute  an  official  Department 
of  the  Army  endorsement  or  approval  of  the  products 
or  services  of  these  organizations. 
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I.  Statement  of  the  Problem 


The  ultimate  goal  of  our  interdiscipl  inary  research  under  this  Army 
contract  is  to  delineate  the  stereoelectronic  and  physicochemical  requisites 
for  oximes  to  be  effective  reactivators  of  acetylcholinesterase  (AChE) 
inhibited  by  organophosphorus  compounds.  This  is  being  done  by  carrying  out 
theoretical  high  quality  three-dimensional  ab-initio  quantum  chemical  cal¬ 
culations  on  the  oximes,  and  from  these  wave  functions  generating  three- 
dimensional  electrostatic  molecular  potential  contour  (EMPC)  maps  around 
these  oximes  as  well  as  experimental  determinations  of  their 
1 ipophil icities. 

The  major  objectives  of  our  research  as  described  in  the  scope  of  work 
of  this  contract  are  as  follows: 

a.  To  carry  out  high  quality  three  dimensional  quantum  chemical  cal¬ 
culations  and  generate  electrostatic  molecular  potential  contour 
maps  of  the  effective  oximes  2-PAM*,  HI-6*  and  at  least  one  other 
oxime  (and  for  various  other  effective  oximes  if  time  and  funds 
permit)  to  delineate  the  stereoelectronic  requisites  necessary  for 
effective  action. 

b.  To  carry  out  similar  quantum  chemical  and  EMPC  map  calculations  for 
oximes  reacting  with  models  of  phosphory 1 ated  AChE  and  of  the  oximes 
themselves  reacting  with  models  of  AChE. 

c.  To  formulate  how  the  quantum  chemical  and  EMPC  map  information  would 
be  used  to  set  up  computer  based  structure-activity  prediction 
program  for  the  efficacy  of  oximes. 

d.  To  use  the  quantum  chemical  and  EMPC  map  information  to  set  up  a 
program  for  lead-directed  synthesis  of  new,  more  effective  oximes. 

e.  To  carry  out  careful  experimental  measurements  of  the  pK  ‘ s  of  the 
oximes  (including  temperature  dependence  and  deconvol ut i$g  overlap¬ 
ping  pK ,’s).  Also,  to  carry  out  experimental  measurements  of  the  pH 
dependence  of  the  1 ipophi 1 icities  (drug  distribution  coefficients) 
of  oxime  reactivators  (including  2-PAM,  HI-6)  using  contractor's 
sensitive  novel  micro-electrometr ic  titration  technique 

to  relate  this  to  observed  in  vivo  activity  to  be  carried  out  else¬ 
where  in  accordance  with  the  direction  of  the  COTR ;  further,  to 
use  these  experimental  measurements  as  a  guide  to  accurate  theoreti¬ 
cal  prediction  of  the  lipophilic  drug  distribution  coefficients. 


*  See  Figure  1  of  this  report  for  sketches  of  the  molecular  structure. 
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II.  Background 

EMPC  maps  around  a  variety  of  drugs  for  various  different  therapies 
(such  as  mentioned  above  in  parts  a-d)  have  previously  been  shown  to  relate 
to  the  effectivities  of  such  drugs.  About  seven  years  ago,  by  invitation. 
Dr.  Kaufman  presented  lectures  at  thirteen  different  pharmaceutical  houses 
both  here  and  abroad,  discussing  our  successes  in  correlation  of  quantum 
chemical  indices  and  EMPC  maps  with  pharmacological  activity  and  in  drug 
design  prior  to  syntheses  with  this  combination  of  theoretical  techniques. 

A  number  of  pharmaceutical  companies  have  now  set  up  theoretical  groups 
attempting  to  use  these  techniques  for  drug  design. 

The  following  discussion  indicates  briefly  how  our  calculated  quantum 
chemical  indices  and  EMPC  maps  are  relevant  to  the  question  of  reactivation 
by  oximes  of  AChE  inhibited  by  organophosphorus  compounds. 

Previous  experimental  in-vitro  pharmacology  studies  have  shown  that 
effectivity  of  the  reactivation  in-vitro  can  be  described  by  two  terms: 
affinity,  K  ,  (due  to  complex  formation  of  the  reactivator  with  the  phos- 
phorylated  AChE)  and  chemical  reactivity,  k  .  The  combination  of  the  two 
terms  is^a  measure  which  has  been  shown  to  correlate  with  the  effectivity 
in-vivo. 

This  two-component  mechanism  of  reactivation  of  inhibited  AChE  has 
conceptual  similarities  to  the  first  two  steps  in  the  mechanism  of  inhibi¬ 
tion  of  AChE  by  phosphorus  compounds,  namely,  formation  of  the 
organophosphorus- AChE  complex  and  then  a  chemical  step  which  involves  break¬ 
ing  a  P-X  bond  and  phosphorylation  of  the  AChE.  We  have  previously  shown 
that  the  first  step  in  inhibition  of  AChE  (the  complex  formation)  is  related 
to  the  EMPC  maps  around  the  organophosphorus  compounds.  We  then  showed 
that  the  second  step  in  the  inhibition  of  AChE  (phosphorylation  of  the  AChE) 
is  related  to  a  theoretical  quantum  chemical  index  of  the'P-X  bond. 

Thus,  by  analogy,  for  the  present  problem  of  reactivation  by  oximes  of 
inhibited  AChE,  we  likewise  expect  the  first  step  of  complex  formation  of 
the  oxime  reactivator  to  the  phosphorylated  AChE  to  be  related  to  the  EMPC 
map  around  the  oximes.  The  second  step  of  chemical  reactivity  is  expected 
to  be  related  to  calculated  quantum  chemical  indices  and  to  the  EMPC  maps. 
[The  oximates--deprotonated  oximes--have  been  suggested  by  experimentalists 
probably  to  be  the  relevant  species.] 

There  is  Army  interest  in  developing  oxime  reactivators  which  will  be 
be  more  effective  in  the  central  nervous  system  (CNS)  than  presently  used 
oximes.  The  more  lipophilic  a  molecules  is  (i.e.,  the  greater  its  propen¬ 
sity  to  partition  from  aqueous  to  lipid  phases)  the  more  effective  it  is  in 
the  CNS  (under  most  usual  physiological  conditions). 
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III.  Approach  to  the  Problem  -  Theoretical 

Thus,  the  first  step  of  our  theoretical  research  has  been  systemati¬ 
cally  to  carry  out  ab-initio  quantum  'hemical  calculations  around  prototype 
oximes  and  from  wave  functions  to  generate  the  EMPC  maps  and  other  theoreti¬ 
cal  indices. 

A.  Pyridinium  oximes  and  bis-pyridinium  oximes 

Due  to  a  lack  of  experimental  molecular  crystal  structures  for  almost 
all  of  these  oximes,  this  problem  of  ab-initio  quantum  chemical  calculations 
of  these  oximes  is  very^time  consuming.  There  was  only  a  very  old  crystal 
structure  for  syn  2-PAM  ,  which  did  not  place  the  H  atoms.  There^had  been  , 
more  recent  experimental  crystal  structure  determinations  for  syn3  and  antib 
4-PA  (4-pyridine  aldoxime).  Thus,  in  order  to  get  appropriate  starting 
structures  for  syn  and  anti  2-PAM,  we  first  carried  out  quantum  chemical 
calculations  on  syn  and  anti  2-PA  (putting  the  aldoxime  side  chain  either 
syn  or  anti  in  the  2-position  of  pyridine).  We  then  carried  out  quantum 
chemical  calculations  on  syn  and  anti  2-PA  as  a  function  of  geometry  (CNOH 
dihedral  angle).  This  gives  us  a  calculated  conformational  energy  profile. 
All  possible  conformations  within  a  certain  number  of  kcal  (on  the  order  of 
12  kcal/mole)  from  the  absolute  energy  minimum  may  be  possible  to  achieve  in 
solution  or  under  the  influence  of  an  active  site,  and  thus  may  be  sig¬ 
nificant  pharmacologically.  We  also  used  differencing  techniques  to  project 
a  minimum  which  we  also  checked  by  calculations.  The  equilibrium  geometry 
of  a  species  is  where  the  total  energy  is  at  a  minimum.  Our  results  of 
quantum  chemical  calculated  total  energy  vs.  dihedral  angle  for  syn  and  anti 
2-PA  showed  that  there  were  multiple  maxima  and  minima  in  these  compounds-- 
thus,  many  points  must  be  computed.  [For  these  ab- initio  quantum  chemical 
calculations,  we  are  using  our  own  fast  ab-initio  progr amsRwh i ch  also 
incorporate  optional  desirable  computational  strategies:'’  ’ 

M0DP0T  -  effective  core  model  potentials  which  allow 
calculation  of  the  valence  electrons  only 
explicitly,  but  accurately; 

VRDDO  -  (variable  retention  of  diatomic  differential 

overlap),  an  efficient  charge  conserving  integral 
prescreening  evaluation  which  tells  reliably 
whether  an  integral  will  be  large  enough  to 
be  significant  or  so  small  that  it  should  not 
be  calculated 

which  cut  down  the  computer  time  by  at  least  an  order  of  magnitude  even  for 
a  compound  with  one  aromatic  ring  and  a  side  group  like  2-PA  and  con¬ 
siderably  more  for  larger  molecules. 

MERGE  -  to  save  invariant  skeletal  integrals  and  reuse  them. 
(This  MERGE  technique  is  especially  effective  when 
geometry  changes  and  also  saves  orders  of  magnitude 
of  computer  time  when  series  of  related  compounds 
are  being  studied) . 
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From  these  calculated  ab-initio  MODPOT/VRDDO/MERGE  wave  functions,  we 
calculate  the  gross  atomic  populations  on  the  atoms  (the  atomic  charge 
distribution  in  the  molecule)  and  other  pertinent  theoretical  indices.] 

Thus,  to  carry  out  the  quantum  chemical  calculations  on  2-PAM,  we 
meshed  together  pieces  of  experimental  crystal  structures  for  various  frag¬ 
ments  and  then  carried  out  the  quantum  chemical  calculations  as  a  function 
of  varying  the  geometry  (CNOH  dihedral  angle)  of  both  the  syn  and  anti  2- 
PAM.  For  an  initial  geometry  for  quantum  chemical  calculations  on  syn  and 
anti  2-PAM,  we  then  ansatzed  a  CPU  group  onto  the  pyridine  ring  of  syn  or 
anti  2-PA  (using  the  exper imenta r crystal  of  a  methyl  pyridinium  as  a 
guide) . 

We  carried  out  the  quantum  chemical  calculations  on  syn  and  anti  2-PAM 
as  a  function  of  geometry  (CNOH  dihedral  angle).  Again,  our  results  of 
total  calculated  energy  vs.  dihedral  angle  showed  that  syn  and  anti  2-PAM 
also  have  multiple  maxima  and  minima  (Figure  2). 

(During  the  course  of  our  investigations,  a  new  crystal  structure 
determination  of  2-PAM  was  reported.  We  will  be  running  more  ab-initio 
MODPOT/VRDDO/MERGE  calculations  on  2-PAM  using  the  new  crystal  structure  as 
a  starting  point.  However,  the  crystal  structure  geometries  are  merely  the 
starting  geometries  for  our  quantum  chemical  calculations.  The  calculation 
are  then  run  as  a  function  of  various  geometry  variations.) 

We  carried  out  similar  quantum  chemical  calculations  for  syn  and  anti 
4-PA  and  syn  and  anti  4-PAM  as  a  function  of  geometry. 

For  the  starting  geometries  of  syn  and  anti, 4-PA,  we  used  the  ex¬ 
perimentally  measured  x-ray  crystal  s tructuresJ  ’  and  then  also  calculated 
the  total  energies  as  a  function  of  CNOH  dihedral  angle.  Again,  there  were 
multiple  maxima  and  minima  in  the  energy  function  of  geometry.  We  carried 
out  the  ab-initio  MODPOT/VRDDO/MERGE  calculations  for  rotamers  of  syn  and 
anti  4-PAM.  For  stajrtgng  geometries  we  used  the  x-ray  crystal  goemetries 
of  syn  and  anti  4-PA3’  plus  a  methyl  group  on  the  pyridine  nitrogen  as  on 
the  methylpyridinium  ion.  Frcm  the  plot  of  total  energy  vs.  CNOH  dihedral 
angle  for  syn  and  anti  4-PAM  (Figure  3)  it  is  apparent  that  again  there  are 
multiple  maxima  and  minima. 

This  aspect  of  research  gave  us  knowledge  about  the  conformational 
complexities  of  mono  pyridinealdoxime  molecules.  As  we  calculated  the 
quantum  chemical  wave  functions  of  those  oximes,  we  generated  electrostatic 
molecular  potential  contour  (EMPC)  maps  in  three  dimensions  around  the 
lowest  energy  conformation  of  each  molecule.  From  these  EMPC  maps  we  gener 
ated  three  dimensional  isopotential  EMPC  maps  around  each  molecule. 
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The  EMPC  maps  are  calculated  point  by  point  for  planes  through  and 
around  the  molecules.  We  customarily  calculate  the  potentials  for  1000 
points  per  plane  and  from  15  to  20  planes  per  molecule  for  molecules  the 
size  of  2-PAM.  From  these  points  we  calculate  three  dimensional  isopoten¬ 
tial  EMPC  maps  for  the  molecules.  All  the  three  dimensional  EMPC  points  and 
isopotential  EMPC  maps  are  stored  in  the  computer  for  future  matching  of  one 
oxime  against  another. 

[Color  copies  of  these  three  dimensional  isopotential  EMPC  maps  were 
included  with  our  Quarterly  Progress  Reports  (QPR's).]  Black  and  white 
copies  of  these  three  dimensional  isopotential  EMPC  maps  are  included  in 
this  report  (Figure  4,  syn  2-PAM;  Figure  5,  anti  2-PAM;  Figure  6,  syn  4-PAM; 
Figure  7 ,  anti  4-PAM) . 

It  had  been  suggested  by  the  developers  of  the  pyridine  aldoximes  that 
the  deprotonated  anions  (the  oximates)  of  these  compounds  were  the  effective 
species  for  the  reactivation  of  the  phosphory 1 ated  (inhibited)  AChE. 

Hence  we  initiated  quantum  chemical  studies  of  the  deprotonated 
(zwitterionic)  forms  of  2-PAM  (the  oximate  species).  Crystalline  2-PAM  had 
been  found  to  be  syn.‘ 

However,  it  had  been  suggested  by  Gray  that. in  the  deprotonated 
species,  the  syn-anti  distinction  is  destroyed.1 

Convergence  is  strongly  dependent  on  the  molecular  geometry.  We  had  to 
carry  out  quantum  chemical  calculations  at  a  large  number  of  geometry  points 
for  deprotonated  2-PAM. 

We  had  no  difficulty  obtaining  a  converged  structure  for  the 
deprotonated  syn  2-PAM  anion  in  the  anti-configuration.  No  converged  syn- 
anions  have  been  obtained  to  date.  In  an  effort  to  get  a  stable  structure 
for  the  2-syn  anion,  we  have  let  the  syn  structure  relax  in  a  number  of 
ways : 

From  the  converged  wave  functions,  deprotonated  (zwitterionic)  anti  2- 
PAM  we  generated  the  EMPC  maps.  [A  color  copy  of  the  three  dimensional 
isopotential  EMPC  map  around  deprotonated  2-PAM  was  included  in  our  QPR  2.] 

A  black  and  white  copy  of  this  three-dimensional  isopotential  EMPC  map  is 
included  in  this  report  (Figure  8). 

We  are  continuing  to  investigate  the  quantum  chemical  calculations  on 
deprotonated  2-PAM  with  other  geometry  variations  (bond  lengths,  bond 
angles,  dihedral  angles,  etc.). 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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We  also  initiated  quantum  chemical  studies  of  the  deprotonated 
(zwitterionic)  forms  of  both  syn  and  anti  4-PAM.  We  have  let  the  structure 
relax  in  a  number  of  ways.  No  converged  anions  of  either  of  the  4-PAM' s 
have  been  obtained  to  date. 

It  is  apparent  from  our  quantum  chemical  studies  on  the  deprotonated  2- 
PAM  anions  (where  after  many  attempts  at  different  geometries  we  were  able 
finally  to  get  a  converged  SCF  wave  function  for  the  deprotonated  anti  2-PAM 
anion)  and  our  quantum  chemical  studies  on  the  deprotonated  4-PAM  anions  in 
the  third  quarter  that  the  convergence  of  these  species  is  very  strongly 
dependent  on  the  geometry  chosen.  A  larger  number  of  different  geometries 
must  be  investigated  by  explicit  quantum  chemical  calculations.  We  are 
continuing  to  investigate  this  problem. 
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While  we  were  continuing  our  quantum  chemical  and  EMPC  studies  on  2-PAM 
and  4-PAM  (and  their  deprotonated  anions),  we  initiated  studies  on  HI-6*,  a 
bis-pyridin ium  monooxime  (which  is  also  an  effective  oxime  reactivator  of 
phosphory lated  AChE)  to  gain  insight  into  the  quantum  chemical  indices  and 
EMPC  maps  of  bis-pyridinium  oximes. 


HI-6 


0^/NH2 


OH 


For  a  starting  geometry  for  the  left  ring  system,  we  used  the  geometry  of  2- 
PAM  as  optimized  in  our  previous  studies;  for  the  right  ring  system,  we.used 
the  geometry  from  the  experimental  crystal  structure  of  nicotinic  acid, 
and  we  used  the  ether  bridge  from  the  crystal  structure  of  dimethyl  ether. 

We  are  investigating  various  rotamers  with  relation  to  the  COCN  dihedral 
angle. 

We  generated  three  dimensional  isopotential  contour  EMPC  maps  around 
one  of  the  rotamers  of  HI-6  (a  bis-pyridinium  monooxime).  (A  color  copy  of 
the  isopotential  EMPC  map  around  HI-6  was  included  on  our  QPR  3.  Black  and 
white  isopotential  EMPC  maps  for  HI-6  are  shown  in  this  annual  report, 

F igure  9. ) 

We  will  be  comparing  these  quantum  chemical  indices  and  isopotential 
contour  EMPC  maps  of  the  pertinent  species  of  HI-6  with  those  of  2-PAM  to 
compare  similarities  and/or  differences.  (During  the  course  of  our 
theoretic^  investigations,  a  crystal  structure  determination  of  HI-6  was 
reported.  We  will  be  carrying  out  our  theoretical  investigations  of  HI-6 
also  using  this  crystal  structure.  However,  we  use  the  experimental  crystal 
structure,  if  available,  only  as  a  starting  point.  We  then  carry  out  the 
quantum  chemical  calculations  as  a  function  of  geometry  variations  to  estab¬ 
lish  the  conformational  energy  profiles  of  the  molecules.) 

We  are  continuing  our  quantum  chemical  and  EMPC  map  calculations  on 
these  pyridinium  oximes. 


★ 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 


B.  Imidazol  ium  oximes  BHH  compounds 

In  June  of  1984,  Dr.  Lily  Tang,  the  Army  Contract  Technical  Monitor  on 
this  project  called  up  and  requested  by  phone  that  we  see  if  our  techniques 
could  shed  insight  on  the  relative  "protective  effects  in  animals  toward 
organic  phosphate  poisoning"  [i.e.  against  inhibition  of  acetylcholines¬ 
terase  (AChE)  by  organic  phosphates]  of  certain  of  the  quaternary 
methyl imidazol ium  oxime  derivatives  being  synthesized  by  Dr.  Scott  Bedford 
at  Stanford  Research  Institute1  (mentioned  in  Dr.  Tang's  letter  to  us  dated 


June  21,  1984).  The  compounds  listed  in 
when  used  together  with  atropine  were 

Code  numbers  in  Dr.  Tang's  letter 

Dr.  Tang's  letter  as  "protective" 

Code  numbers  as  assigned  to  us  by 

Col 

.  Pick  our  original  Army  COTR 

BK64397 

BHH-0013 

BK64404 

BHH-0014 

BK64413 

BHH-0015 

BK64422 

BHH-0016 

BK64431 

BHH-001 7 

The  relative  protectivities  as  confirmed 

in 

our  subsequent  telephone  conver 

sation  with  Dr.  Tang  were 

Compound 

Dose  (i.m.) 

BK64397 

BHH- 1 3 

1.8  mg/Kg 

BK64404 

BHH- 14 

7.2  mg/Kg 

BK64413 

BHH- 1 5 

16.7  mg/Kg 

BK64422 

BHH - 1 6 

21.9  mg/Kg 

BK64431 

BHH- 1 7 

33.3  mg/Kg 

Dr.  Tang  also  asked  if  our  technique  could  set  up  theoretical 
stereoelectronic  requisites  which  could  help  distinguish  between  BHH  com¬ 
pounds  which  had  effectivity  as  reactivators  of  inhibited  phosphorylated 
AChE  and  those  which  had  "protective"  effects. 

This  was  a  major  change  in  our  research  plans  under  this  Army  contract 
and  also  greatly  increased  the  amount  of  theoretical  research  to  be  done. 

To  be  responsive  to  Dr.  Tang's  request  we  turned  our  attention  to 
considering  this  problem  of  the  quaternary  methyl imidazol ium  oximes.  There 
were  several  aspects  which  had  to  be  considered: 

1.  These  particular  BHH  "protective"  compounds  were  larger  and  had 
more  conformational  flexibility  than  the  oximes  which  were  agreed  upon 
jointly  by  the  Army  Medical  Research  and  Development  Command  and  by  us  in 
the  scope  of  work  of  our  contract  for  quantum  chemical  investigations. 

These  BHH  compounds  also  increased  greatly  the  number  of  oximes  to  be  inves¬ 
tigated  theoretically. 


t^t  l».ji 


I  ft**  *-*  ^  J 


To  calculate  BHH  compounds  necessitated  our  decision  to  devote  our 
immediate  attention  to  expediting  our  changeover  to  the  new  vector  supercom¬ 
puters  for  quantum  chemical  calculations.  In  order  to  transmit  the  large 
data  files  from  the  calculations  back  to  our  labs  at  the  Johns  Hopkins 
Univeristy,  we  are  having  to  establish  and  implement  asynchronous  error 
detecting,  error  correcting  protocols.  Which  machine  will  be  used  for  which 
particular  problems  will  depend  on  the  core  memory  (virtual  memory,  if  any) 
and  disc  space  available  at  any  particular  time.  This  is  important  when 
molecules  start  to  get  large--of  the  size  of  some  of  the  BHH  molecules. 

2.  We  initiated  ab-initio  MODPOT/VRDDO  quantum  chemical  calcula¬ 
tions  and  electrostatic  molecular  potential  contour  (EMPC)  map  generation  on 
the  quaternary  methyl imidazol ium  oxime  derivatives. 

Since  these  quantum  chemical  calculations  must  be  carried  out  for 
several  such  BHH  derivatives  and  we  will  be  using  the  MERGE  technique,  we 
are  carrying  out  these  quantum  chemical  calculations  systematically  starting 
with  the  parent  compounds.  To  date  we  have  carried  out  the  theoretical 
calculations  for  imidazole,  imidazol ium  (protonated  imidazole)  and  BHH-I. 

BHH-1  has  been  reported  by  Dr.  Bedford  of  SRI  to  be  a  good  reactivator 
in-vitro  against  eel  AChE  inhibited  by  EPMP  (ethy Iparanitropheny 1  methyl 
phosphonate)  or  by  GD  (soman). 

For  imidazoJe  itself  there  was  a  crystal  structure  geometry  to  use  as  a 
starting  point.1  The  quantum  chemical  and  EMPC  map  calculations  on  im¬ 
idazole  also  have  an  intrinsic  interest  with  respect  to  reactivation  of 
phosphorylated  AChE  because  very  early  experimental  reports  indicated  that 
the  inhibited  enzyme  can  be  dephosphorylated  and  thus  reactivated  bye- 
nucleophilic  reagents  such  as  imidazole,  pyridine  and  hydroxyl  ami ne10 ’ 
(reactivation  by  these  reagepts7is,  however,  rather  slow,  so  that  high 
concentrations  are  required10’  ). 

We  also  carried  out  ab-initio  MODPOT/VRDDO  calculations  for  imidazolium 
(protonated  imidazole)  and  generated  the  EMPC  maps.  There  was  also  a  crys¬ 
tal  structure  for  imidazolium  to  use  as  a  starting  geometry.1 

There  is  no  experimental  crystal  structure  for  BHH-1.  Thus  for  a 
starting  geometry  we  constructed  a  hybrid  structure  meshing  together. ^n 
N-methyl  imidazole  fragment  (for  which  there  was  a  crystal  structure1  j  with 
a  syn  aldoxime  s4de  chain  (from  the  crystal  structure  of  syn  4- 
pyr idinealdoxime0)  and  with  a  -CH--0-CH,  side  chain  on  the  other  N  (using 
the  crystal  structure  parameters  from  dimethyl  ether). 

We  carried  out  ab-initio  MODPOT/VRDDO  calculations  for  BHH-1. 


[See  Figure  next  page] 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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BHH-1  CH, 

i  j 


ch2och3 


We  are  contin 
dihedral  angles, 
around  one  of  the 


uing  the  calculations  as  a  funct 
We  generated  a  three  dimensional 
rotamers  of  BHH-1  (Figure  10). 


on  of  geometry  and 
isopotential  EMPC  map 


We  plan  to  compare  the  EMPC  maps  for  2-PAM  and  for  BHH-1  (and  their 
deprotonated  species)  to  identify  similarities  [or  possibly  differences 
(while  since  both  are  good  react  i  vators ,  these  compounds  have  somewhat 
different  reactivating  properties)]  between  the  stereoelectronic  charac¬ 
teristics  of  these  two  molecules.  We  plan  to  continue  carrying  out 
calculations  for  these  BHH  quaternary  methyl imidazol ium  oxime  compounds  to 
the  extent  our  Army  contract  resources  permit. 


3.  Possible  mechanisms  contributing  to  protective  actions  in-vivo 
of  certain  of  these  BHH  oximes 

At  our  request,  Dr.  Tang  got  Dr.  Bedford  to  send  us  his  technical 
reports  on  the  BHH  quaternary  methyl imidazol ium  oximes.  We  received  from 
him  a  preprint  manuscript.  We  examined  this  manuscript  carefully.  While 
this  manuscript  contained  data  on  the  reactivation  in-vivo  of  phosphory lated 
AChE  by  these  BHH  compounds  there  was  no  mention  of  the  protective  in-vivo 
animal  experiments  to  which  Or.  Tang's  letter  had  referred.  However,  upon 
reading  that  preprint  manuscript,  it  appeared  that  the  majority  of  BHH 
compounds  which  were  effective  reactivators  were  different  from  those  which 
were  "protective" . 

Since  there  was  no  direct  experimental  data  in  Dr.  Bedford's  report 
which  could  shed  insight  into  the  mechanistic  aspects  of  the  protective 
effect  of  the  BHH  oximes,  it  was  necessary  to  find  out  from  previous 
reported  experimental  results  what  mechanisms  could  be  contributing  to 
"protective  action"  of  oximes  and  consider  these  for  the  BHH  oximes. 


*  See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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We  then  ascertained,  from  examining  prior  experimental  research 
results,  that  compounds  which  had  prophylactic  activity  against  inhibition 
of  AChE  by  organophosphorus  compounds  could  do  so  by  several  very  different 
mechanisms.  (Also  prophyl active  activity  was  dependent  on  the  chemical 
nature  of  the  reactivator  and  of  the  particular  inhibitor.) 

a.  Compounds  protective  against  organophosphorus  inhibition  of 
AChE  can  react  with  the  AChE  active  site  (either  with  the  catalytic  es- 
teratic  site  and/or  with  the  neighboring  anionic  site)  and/or  with  the 
active  site  and  with  the  peripheral  anionic  site^gc  even  just  with  the 
peripheral  anionic  site  (if  sufficiently  bulky).  Moreover,  interaction  at 
the  peripheral  AChE  anionic  site  can  change  the  conformation  of  the  active 
AChE  site. 


b.  Compounds  protective  against  organophosphorus  inhibition  of 
AChE  can  react  directly  in-vivo  with  the  organophosphorus  AChE  inhibitor. 

HS-6  has  been  reported  to  react  with  soman  under  physiological 
conditions  by  accelerating  the  hydrolysis  of  soman  before  it  is 
able  to  react  with  AChE  .  HS-6  shows  protective  effects  in-vivo. 

c.  Certain  oxime  reactivators ,  such  as  2-PAM,  do  not  prevent  an 
increase  in  brain  AChE  after  exposure  to  organophosphorus  AChE  inhibitors 
(2-PAM  itself  prophy lactical ly  actually  causes  a  rise  in  brain  AChE)  but  are 
protective  against  symptoms.  This  has  been  suggested  as  being  due  to  the 
cholinergic  receptors  becoming  tolerant  to  increased  levels  of  ACh. 

io  date  we  have  not  yet  received  any  other  information  from  Dr.  Tang  or 
Dr.  Bedford  on  experimental  results  which  could  clarify  these  points  con¬ 
cerning  the  mechanism  of  protective  effects  of  the  BHH  compound. 

From  our  analysis  there  are  several  specific  other  experiments  which 
might  shed  light  on  some  of  these  points.  Thus,  in  the  absence  of  ex¬ 
perimental  mechanistic  data  we  plan  to  look  for  differences  in  our 
calculated  quantum  chemical  EMPC  map  results  between  the  BHH  quaternary 
methyl imidazol ium  oximes  which  are  effective  as  reactivators  and  those  which 
are  effective  as  "protective"  agents. 

C.  Possible  steps  in  mechanism  of  reactivation  of  phosphory 1 ated  AChE 

One  of  the  objectives  of  our  original  scope  of  work  was: 

b.  To  carry  out  similar  quantum  chemical  and  EMPC  map  calculations 
for  oximes  reacting  with  models  of  phosphory 1 ated  AChE  and  of 
the  oximes  themselves  reacting  with  models  of  AChE. 

Thus,  we  began  to  examine  what  might  be  appropriate  approaches  for 
understanding  the  mechanism  of  reactivation  of  phosphory 1 ated  AChE  to  derive 
appropriate  quantum  chemical  and  other  indices. 

It  had  been  shown  that  effectivity  of  the  reactivation  in-vitro  can  be 
described  by  two  terms:  affinity,  K  ,  (due  to  complex  formation  of  the 
reactivator  with  the  phosphory 1 ated  AChE)  and  chemical  reactivity,  k  .  The 
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combination  of  the  two  terms  is  a  measure  which  has  been  shown  to  correlate 
with  the  effectivity  in-vivo.1 

The  Bedford  manuscript  did  contain  K  and  kr  results  for  the  reactiva¬ 
tion  of  EPMP-inhibited  eel  AChE  by  the  BHfi  compounds,  but  the  aging  of 
soman- inhibited  AChE  was  sufficiently  rapid  to  preclude  reliable  measurement 
of  K  and  k  ,  so  the  maximum  reactivation  of  AChE  was  used  by  them  as  a 
measure  of  reactivating  efficacy. 

We  plan  to  examine  reactivation  of  phosphoryl ated  AChE  from  the  above 
point  of  view  of  affinity  and  of  chemical  reactivity  since  our  prior  ex¬ 
perience  on  quantum  chemical  calculations  and  EMPC  map  generation  for  the 
inhibition  of  AChE  by  organophosphorus  compounds  indicated  that  each  of  the 
three  separate  steps  in  the  inhibition  mechanism  of  AChE  was  governed  by  a 
different  theoretical  index. 

As  a  measure  of  the  affinity  factor,  we  are  examining  the  size  and 
shape  of  the  positive  and  negative  EMPC  contours  around  the  deprotonated  2- 
PAM  (color  copies  of  which  were  included  in  our  QPR2;  black  and  white  copies 
are  included  in  this  annual  report,  Figure  8)  with  respect  to  the  inferences 
that  have  been  made  previously  for  the  spatial  relation  between  the  es- 
teratic  and  anionic  subsites  of  the  AChE  active  site. 
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IV.  Results  and  Discussion  of  the  Results  -  Theoretical 

The  three  dimensional  isopotential  contour  EMPC  map  for  the 
deprotonated  2-PAM  showed  vividly  a  most  significant  finding.  [A  color  copy 
of  the  isopotential  EMPC  map  around  deprotonated  2-PAM  was  included  with  our 
Quarterly  Progress  Report  2  (QPR  2).  A  black  and  white  isopotential  EMPC 
map  for  deprotonated  2-PAM  is  shown  in  this  Annual  Repo  t  as  Figure  8.]  The 
shape  and  disposition  of  isopotential  contours  of  this  map  indicates  vividly 
that  although  the  deprotonated  2-PAM  is  an  overall  neutral  species,  [since 
there  is  a  unit  (+)  charge  from  the  quaternization  of  the  nitrogen  and  a 
unit  (-)  charge  from  the  deprotonation  of  the  oxime  H],  the  isopotential 
EMPC  map  indicates  this  species  to  be  a  dipolar  ion.  The  three  dimensional 
isopotential  EMPC  map  of  deprotonated  2-PAM  indicated  clearly  that  there  is 
a  strongly  positive  region  and  a  strongly  negative  region.  This  has  sig¬ 
nificant  implications  for  several  aspects  of  the  use  of  oxime  reactivators. 
First,  the  size  and  shape  and  relative  spatial  relationship  of  the  positive 
and  negative  isopotential  portions  of  these  EMPC  maps  should  have  a  com¬ 
plementary  relation  to  the  spatial  relationship  and  electronic 
characteristics  of  esteratic  and  anionic  subsites  of  the  active  site  of 
AChE.  The  is  particularly  important  in  view  of  the  complete  lack  of 
molecular  crystal  structural  information  on  any  AChE.  Our  isopotential  EMPC 
maps  of  deprotonated  2-PAM  can  help  shed  insight  onto  the  size,  shape  and 
stereoel ectronic  characteristics  of  the  AChE  active  sites--quantities  at 
present  not  yet  obtained,  since  so  far  it  has  apparently  not  been  possible 
to  crystallize  AChE  to  determine  its  crystal  structure. 

Second,  the  isopotential  EMPC  map  of  deprotonated  2-PAM  is  as  unique  as 
a  fingerprint  in  identifying  stereoelectronic  characteristics  that  have  the 
proper  stereoelectronic  complementarity  and  are  favorable  for  effective 
reactivation  of  AChE  inhibited  by  organophosphorus  compounds.  The  cor¬ 
responding  calculated  isopotential  contour  EMPC  maps  of  the  oximates  of 
proposed  hypothetical  new  oxime  reactivators  can  be  computer-compared  to 
those  of  deprotonated  2-PAM  prior  to  the  synthesis  of  the  new  compounds  to 
see  if  these  new  compounds  have  the  proper  stereoelectronic  character istics 
to  be  effective  new  reactivators. 

Third,  our  three-dimens ional  isopotential  EMPC  maps  around  deprotonated 
2-PAM  also  shed  insight  into  another  important  aspect  of  the  properties  of 
2-PAM  and  deprotonated  2-PAM.  2-PAM  itself  is  a  positively  charged  quater¬ 
nary  species.  As  such,  it  is  not  expected  to  be,  nor  is  it,  very 
lipophilic.  2-PAM  has  little  tendency  to  partition  from  aqueous  to  lipid. 

It  had  been  hoped  that  when  2-PAM  was  deprotonated  and  was  then  an 
overall  neutral  species,  that  it  would  have  a  much  greater  tendency  to 
partition  into  lipid. 

Our  experimental  determination  of  the  1 ipophi 1 ic i ty  of  2-PAM  as  a 
function  of  pH  indicated  that  even  when  the  positively  charged  quaternary  2- 
PAM  was  deprotonated,  the  resulting  species  was  still  only  very  slightly 
lipophilic.  This  behavior  is  commensurate  with  our  EMPC  map  results, 
namely,  that  deprotonated  2-PAM  (while  overall  neutral)  has  the  electronic 
characteristics  of  a  strongly  dipolar  ion. 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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Because  2-PAM  would  have  little  tendency  to  partition  through  the 
lipophilic  blood-brain  barrier  under  usual  physiological  conditions,  2-PAM 
is  not  an  effective  reactivator  in  the  CNS  of  AChE  inhibited  by  organophos- 
phorus  compounds. 

Our  EMPC  map  of  deprotonated  2-PAM  indicated  that  while  this  species  is 
neutral  overall,  it  appears  to  be  a  strongly  dipolar  ion  and  this  would  not 
be  expected  to  have  a  significant  tendency  to  partition  into  lipids.  One  of 
the  greatest  values  of  quantum  chemical  calculations  and  EMPC  maps  is  that 
these  can  be  calculated  for  an  hypothetical  molecule  prior  to  to  its 
synthesis.  If  the  calculated  quantum  chemical  indices  and  the  EMPC  maps 
indicate  that  a  particular  new  oxime  (or  other  proposed  reactivator)  has 
favorable  stereoel ectron ic  characteristics,  then  it  may  be  desirable  to 
pursue  its  synthesis  and  testing.  If,  on  the  other  hand,  the  quantum  chemi¬ 
cal  indices  and  EMPC  maps  indicate  that  a  molecule  does  not  seem  to  be 
favorable  (i.e.,  stereoel ectron ic  character i sti cs  of  its  EMPC  maps  indicate 
that  it  will  not  be  compl ementary  to  the  AChE  active  site  or,  for  CNS  reac¬ 
tivation,  that  its  oximate  may  be  a  strongly  dipolar  ion  with  little 
tendency  to  partition  into  the  CNS  under  specified  blood-brain  barrier 
conditions),  then  the  synthesis  and  testing  of  the  molecule  would  be  less 
desirable  as  a  candidate  for  an  effective  CNS  reactivator. 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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V.  Approach  to  the  Problem  -  Experimental 

In  addition  to  our  theoretical  work,  we  have  been  carrying  out  ex¬ 
perimental  determinations  of  the  1 ipophil icities  of  these  oxime 
reactivators.  We  have  carried  out  these  1 ipophi 1 icity  determinations,  first 
with  our  microelectrometric  technique.  It  turns  out  that  neither  these 
quaternary  oximes  nor  their  deprotonated  species  are  very  lipophilic. 
However,  since  there  have  been  reports  that  some  oxime  reactivators  are 
effective  at  low  concentrations,  we  have  set  up  and  are  also  using  a  more 
sensitive  UV  spectrophotometr ic  procedure  to  establish  quantitatively  and 
directly  just  how  much  of  these  compounds  is  partitioning  into  the  lipid 
phase. 

A.  Microelectrometric  titration  technique 

This  microelectrometric  titration  technique  enables  one  to  determine 
not  only  the  dissociation  constant  of  compounds  which  may  have  more  than  one 
pK  (by  deconvolution  the  overlapping  curves)  but  also  the  solubilities  and 
oit-water  (n-octanol /water)  distribution  of  compounds  as  a  function  of  pH. 

The  octanol/water  partition  coefficients  (P)  of  the  unionized  compound 
were  determined  by  microelectrometric  titration  by  measuring  the  apparent 
shift  in  the  pK  value  in  the  presence  of  octanol.  In  the  aqueous  phase, 
the  compound  dissociates  into  its  respective  ions  according  to  its  K&  value: 

Ka  =  [A]  [H+]/[HA] 

Where  the  dissolved  compound  is  titrated  in  the  presence  of  finely 
dispersed  octanol,  the  protonated,  undissociated  compound  distributes  be¬ 
tween  the  octanol  and  water.  The  removal  of  HA  into  the  octanol  results  in 
a  shift  in  the  apparent  pK  ,  by  an  amount  which  depends  on  the  volume  (V^) 
of  octanol  and  the  volume  fv  )  of  water.  From  this  shift,  the  partition 
coefficient  (P)  can  be  calculated: 


P  = 


1] 


The  partition  coefficient  (P)  is  defined  as  the  partitioning  of 
precisely  the  same  species  between  oil  and  water,  i.e. 


Woil 

^^water 


However,  the  total  drug  distribution  coefficient  (D)  at  a  certain  pH, 
which  is  pertinent  to  an  understanding  of  the  defferences  in  pharmacologica 
potency,  is  the  ratio: 

EAH  ^  Opil  -  Woil 

£AH  +  A"  water  CAH  +  A"Wer 

[Customarily,  we  determine  the  1 ipophi 1 icities  of  drugs  at  37° C  (which 
is  physiological  body  temperature)  and  sometimes  at  20° C  which  is  a  fairly 
standard  temperature  at  which  the  1 ipophi 1  icities  of  various  chemicals  have 
been  reported.  At  the  request  of  Dr.  Lily  Tang,  the  Army  Contract  Monitor 
on  this  project,  we  also  initiated  studies  at  30° C  for  these  oximes. 
However,  after  initial  studies  at  30° C  on  several  oximes,  fir  the  oximes 
were  not  lipophilic  at  37° C,  we  did  not  carry  out  further  1 ipophi 1 icity 
determinations  at  30° C.] 


B.  UV  spectrophotometric  technique 


The  octanol /water  distribution  coefficients  can  be  determined  from  the 
electrometr ic  titration  data  if  the  titration  is  carried  out  in  a 
water/octanol  mixture  as  described  above. 

However,  this  electrometric  titration  method  is  not  a  suitable  one  if 
the  distribution  coefficient  is  very  small.  Since  this  has  turned  out  to  be 
the  case  with  many  of  the  oximes  we  have  under  study,  we  developed  the 
following,  more  sensitive,  spectrophotometr ic  approach  using  ultraviolet 
absorption  spectroscopy.  The  concentrations  of  the  oximes  in  both  aqueous 
and  octanol  phases  were  determined  by  UV  absorption.  Frequently,  in  using 
this  technique,  the  concentration  in  the  aqueous  phase  is  measured  and  the 
concentration  in  the  oil  phase  is  determined  by  difference  since  the  total 
amount  of  material  used  is  known.  This  latter  approach  was  deemed  not 
completely  suitable  for  the  oximes  since  the  amount  of  material  in  the 
octanol  phase  is  very  small  and  large  errors  can  result. 

The  concentrations  of  the  oxime  in  both  phases  were  determined  from  UV 
absorption  from  the  equation 

A  =  ztc 

where  A  =  optical  absorbance,  Z  =  molar  absorptivity  of  the  solute,  t  is  the 
path  length  in  cm  and  c  is  the  molar  concentration  of  the  solute. 

For  each  of  the  compounds  studied,  the  absorption  maximum  (a  )  and  z 
were  determined  for  the  different  ionic  species  in  aqueous  solution  at  a  pH 
=  1.0  (0%  oxime  anion),  at  a  pH  =  13  (100%  oxime  anion)  and  at  physiological 
pH  =  7.4.  A  calibration  curve  (absorbance  vs  concentration)  was  also  deter¬ 
mined  for  the  octanol  solution  to  establish  that  Beer's  Law  is  obeyed  in  the 
working  concentration  range.  Finally,  checks  were  made  of  the  hydrolytic 
stability  of  the  oximes  at  3  7  0  C  and  pH  =  7.4.  The  hydrolysis  of  the  im- 
idazolium  compounds  was  found  to  be  negligible  for  times  of  up  to  six  hours. 

The  samples  were  mixed  and  put  into  a  shaker  at  37 °C  for  an  appropriate 
period  of  time.  The  two  phases  were  separated  by  centrifugation  and  the 
amount  of  material  in  each  phase  was  determined  from  the  calibration  curves. 
In  each  case,  good  material  balance  was  obtained  and  all  determinations  were 
carried  out  in  triplicate.  The  distribution  coefficient  (D)  was  determined 
from  the  relation 


where  C,  =  concentration  in  octanol  phase,  V,  =  volume  of  octano1,  C2  = 
concentration  in  buffer  phase  and  V2  =  volume  of  buffer.  The  pat  M* ion 
coefficient  was  calculated  from  the  relation 

P  =  D  [1  ♦  10(pH  '  pKa}] 

where  P  is  partition  coefficient  and  0  is  drug  distribution  coefficient. 

The  procedure  is  time  consuming;  however,  when  small  1 ipophi 1 icities 
are  involved,  it  is  a  method  of  necessity. 
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C.  Oximes  investigated 

To  date  we  have  carried  out  microelectrometric  1  ipoph i 1  ic  i  ty  determina¬ 
tions  on: 


2-PAM 

BHH-1 

TMB-4 

BHH-3 

HI-6 

BHH-4 

T0X0G0N  IN 

BHH-5 

BHH-3 

and  UV  spectrophotometr ic  1 ipoph i 1 ic i ty  determinations  on: 

2-PAM 

BHH-13 

BHH-14 

BHH-15. 

The  detailed  results  have  been  reported  on  our  QPR's  and  are  included 
in  the  next  section  of  this  report. 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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VI.  Results  and  Discussion  -  Experimental 
A.  Lipophi 1 icity  determinations 

1.  2 -PAM 

a.  Microel ectrometr ic  technique 

The  ionization  constants  and  the  octanol/water  partition  coefficients 
were  determined  for  2-pyridine  aldoxime  choloride  (2-PAM)  at  20°C,  30°,  and 
37°  C . 


The  following  observations  were  made: 

The  titration  curves  (after  the  subtraction  of  blank),  show  only  one 
titratable  group,  i.e.,  the  oxime  moiety. 

The  aqueous  solution  of  2-PAM  in  water  is  colorless.  During  the  titra¬ 
tion  of  the  sample  with  NaOH,  the  color  of  the  solution  turns  yellow.  The 
sample  with  HC1  follows  the  same  course  in  reverse,  indicating  that  the 
formation  of  the  yellow  colored  derivative  is  an  entirely  reversible 
process.  At  the  end  of  titration  (pH  =  5.5),  the  solution  is  colorless 
again. 

Below  is  the  summary  of  the  results.  The  pK  values  represent  the 
average  of  at  least  4  determinations.  a 


TEMP 

CONDITIONS 

pKa 

20°  C 

Aqueous 

7.904  ± 

0.01 

n-octanol /water 

7.895  ± 

0.01 

30°  C 

Aqueous 

7.785  ± 

0.01 

n-octanol /water 

7.773  ± 

0.01 

37°  C 

Aqueous 

7.675  ± 

0.01 

n-oct  -.nol/water 

7.674  ± 

0.01 

Since  the  difference  in  the  pK  '  s  determined  in  aqueous  or  n- 
octanol /water  solution  is  within  experimental  error,  we  conclude  from  our 
microelectrometric  determinat ions  that  to  within  experimental  error  there  is 
virtually  no  partitioning  of  any  of  the  species  of  2-PAM  from  water  to  oil 
at  any  pH  in  the  range  from  2CP  C  to  37°  C. 

Hence,  neither  2-PAM  (which  is  a  positively  charged  ion  and  does  not 
appear  from  these  microelectrometric  experiments  to  partition)  nor  the 
deprotonated  anion  of  2-PAM  (which  while  formally  neutral,  still  behaves  as 
a  dipolar  ion  and  does  not  partition  from  water  to  octanol)  appear  to  have 
much  tendency  to  pass  through  any  lipophilic  barriers  including  the 
blood/brain  barrier. 

A  brief  survey  of  the  literature  revealed  the  following  pK  values 
assigned  for  2-PAM.  a 


*  See  Figure  I  of  this  report  for  sketches  of  the  molecular  structures. 
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Ginsburg  and  Wilscm  deterajjned  a  pK  at  25  C  by  measuring  the  pH  of  a 
half  neutralized  solution;  Mason^3  reportsaa  value  of^S.OO  ±  0.01  at  20 °C 
measured  by  potentiometr ic  titration;  Hagedorn  et  al.  report  a  value  of 
7.68  i  0.03  by  potentiometr ic  titration  (no  temperature  stated). 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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Since  some  authors  indicate  that  the  reactivation  efficiency  of  ACh 
is  related  to  the  concentration  of  the  oxime  anion  at  physiological 
pH,  we  calculated  the  function  of  the  anion  at  the  range  7.10  -  7.80 
from  the  relationship: 

_  antilog  (pKa  -  pH) 

a  '  -  =  (for  2 -PAM  at  37 ?) 

antilog  (pKa  -  pH)+l 


pH 

%  oxime 

6.80 

11.8 

7.10 

21.0 

7.20 

25.1 

7.35 

32.1 

7.40 

34.7 

7.45 

37.3 

7.50 

40.1 

7.60 

45.7 

7.675 

50.0 

7.70 

51.4 

This  night  provide  especially  useful  in  the  case  of  dioximes, 
where  the  true  proportion  of  the  oxime  anion  to  molecules  at 
pnys iological  pH  can  be  calculated  knowing  the  two  ionization  constants. 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 


b.  UV  spectrophotometr ic  technique 

2-PAM  was  one  of  the  oximes  for  which  the  electrometric  measurements 
indicated  a  very  small  distribution  coefficient  between  octanol  and  water. 

We  began  by  first  investigating  the  UV  spectra  of  2-PAM  in  HC1,  NaOH 
and  buffer  at  pH  =  7.4,  in  order  to  establish  the  molar  extinction,  at  A  = 
293  nm  and  from  there  a  calibration  curve  for  different  concentrations. 

Z  =  1.04  x  104  at  a  =  293  nm 

which  means  that^the  lower  limit  of  the  detectable  amount  of  2-PAM  in  oc¬ 
tanol  is  2  x  10  °  g/1. 

c; 

We  determined  the  solubility  of  2-PAM  in  octanol,  which  is  2  x  10  g/1 

(1.2  x  10'4M). 


The  following  characteristic  constants  were  determined  for  2-PAM  from 
the  optical  measurements: 


SOLUTION 

pH 

%  OXIME  ANION 

A  (nm) 

max 1  1 

Z 

0 . IN  HC1 

1 

0% 

293 

1.219 

X 

o 

0 . IN  NaOH 

13 

100% 

336 

1.310 

x  104 

Phosphate 

Buffer 

7.4 

34.7% 

295 

1.064 

^r 

o 

X 

30 


For  all  the  following  microelectrometr ic  studies,  the  determinations 
were  carried  out  at  37°C  by  titration  in  a  10%  oil/water  solution.  Since 
the  experimentally  measured  pK  values  at  37 °C  were  the  same  to  within 
experimental  error  in  aqueous  solution  or  in  10%  oil/water  solution,  the 
1 ipophi 1 icities  of  these  compounds  were  all  very  small  at  37  C  and  were 
expected  to  be  very  small  at  any  lower  temperature.  Hence  no  further 
1 ipophi 1 icity  determinations  were  made  on  such  compounds  at  lower 
temperatures. 


2.  HI-6 


Microelectrometric  technique 
The  results  are  the  following: 


Temperature 


Conditions 


37  °C 

Aqueous 

7.082 

1  0.01 

37  °C 

10%  oil/water 

7.081 

1  0.01 

y* 

V* 

20  °C 

Aqueous 

7.302 

1  0.014 

J 

The  following  observations  were  made: 


The  titration  curve  (after  the  subtraction  of  the  blank  curve)  shows 
only  one  titrable  group,  i.e.  the  oxime  moiety. 


-  The  aqueous  solution  of  HI-6  (cone.  8  x  10  M)  is  pale  yellow  even 
at  pH  =  5.  The  yellow  color  intensifies  gradually  with  the  increase 
of  pH. 


At  37  C  the  compound  starts  to  hydrolyze  at  pH  values  above  9.75. 
The  retitration  curve  (with  HC1)  does  not  retrace  the  original  titra¬ 
tion  curve  and  the  final  solution  has  a  deep  yellow  color  even  <^t  pH 
=  2.50.  The  compound  hydrolyzes  40-45%  in  ca  15  minutes  at  37  C  at 
pH  >  10.00. 


At  20 °C,  the  hydrolysis  process  (starting  at  higher  pH)  does  not 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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interfere  with  the  pKa  determination  of  the  oxime  group. 


-  The  pKa  of  the  amide  group  lies  well  outside  the  range  of  potentio- 
metric  determination  and  physiological  pH. 

-  The  value  reported  in  literature  is  pKa  =  7.28,  no  temperature  sta¬ 
ted^. 

Since  it  has  been  postulated  that  the  anion  is  the  effective  reactivating 
species,  we  calculated  the  fraction  of  oxime  anion  in  the  physiological  pH 
range  from  pH  =  6.80  to  pH  =  7.70,  from  the  relationship 

anti  log  (pKa  -  pH) 
a  '  anti  log  (pKa  -  pH)  +  1 


:  Oxime  amor.  I 


20°  pKa  =  7.30 


37°  pKa  -  7.08 


24.0 

34.4 

38.7 

51.2 

44.3 

56.9 

50.0 

62.4 

As  one  can  see  from  the  table,  the  temperature  makes  a  difference.  At 
pH  =  7.40,  the  normal  physiological  pH,  the  percentage  of  the  oxime  ion 
increases  217  from  20°C  to  37°C. 

-HI-6  did  not  partition  between  oil/water  at  37’,  thus  no  determinations 
were  made  at  20°C  in  oil/water  solution  since  it  is  not  expected  that 
the  compound  will  partition  into  the  oil  phase  at  the  lower  temperature. 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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3.  TMB-4 

Microel ectrometric  technique 


The  determinations  were  made  at  37°C  and  20"C.  The  titration  curve,  af 
ter  subtraction  of  the  blank  curve,  shows  two  titrable  groups,  i.e.  the  two 
oxime  moieties. 

The  results  are  as  follows: 

Temperature  Conditions  pKai  pKa2 


37°C 

37°C 

20°C 


Aqueous 
100  oil /water 
Aqueous 


7.61  i  0.02 
7.60  i  0.02 
7.82  r  0.03 


8.33  -  0.02 
8.35  t  0.02 
8.59  i  0.03 


23 

Our  previously  proven  technique  of  grapnically  deconvol uti ng  overlap¬ 
ping  pKa's  was  applied  for  the  determination  of  the  two  pKa's. 

From  graphical  deconvolution  the  two  pKa's  are  only  0.72  pH  units  apart. 
That  means  that  the  ionization  of  the  second  oxime  group  starts  even  before 
half  of  the  first  oxime  group  is  ionized. 


The  literature  value  for  TMB-4  is  pK  =  7.78  -  0,01,  pfC  =  8.61  •  0.03; 

27  28  31  aZ 

no  temperature  given.  ’ 

The  fraction  of  the  oxime  anion  in  the  physiological  significant  range 
is  the  sum  of  the  mono  aldoxime  and  bispyridiniumaldoxime  anions. 

The  compound  does  not  hydrolyze  readily  at  pH  up  to  11.50  in  15  min., 
even  at  37 'C  --  the  retitration  curve  with  HC1  retraces  the  original  titration 
curve. 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures 
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Since  the  two  pK  values  lie  again  quite  close  to  one  another-- 1  ess 
than  one  pK  unit  apar?--our  technique  of  graphical  deconvolution  of  overlap¬ 
ping  pK  's  may  be  less  accurate  than  would  be  desirable. 

a 

Therefore,  we  set  up  a  computer  progra^gfor  calculating  pK^  values 
based  on  an  equation  developed  by  Speakman. 

However,  for  the  separate  pK  's,  our  graphical  results  and  the  results 
from  the  computer  program  were  virtually  identical. 


♦ 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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4.  TOXOG0NIN 

The  determinations  were  made  at  37°C  and  20°C.  The  results  are  the 
following. 


Temperature 

Condi tions 

pKai 

pKat 

37°  C 

Aqueous 

7.30  ±  0.02 

8.14  ±  0.02 

10%  oil  /  water 

7.30  -  0.04 

8. 17  r  0.01 

20°  C 

Aqueous 

7.55  i  0.02 

8.32  t  0.04 

All  the  observations  made  for  TMB-4  apply  to  T0X0G0NIM  also.  We 
used  the  graphic  method  for  deconvoluting  the  two  overlapping  pKa's. 

The  compound  does  not  hydrolyze  at  37°  in  about  15  minutes  at  pH  up  to 
11.50.  The  retitration  curve  with  HC1  retraces  the  original  titration 
curve. 


The  literature  values  are  7.54  :  0.01  and  8.24  r  0.2  -  no  tempera¬ 
ture  gi ven .  29 

Hagedorn25! ists  a  value  of  7.80,  which  she  derived  from  pKa}  =  7.54 
and  pKa„  =  8.12,  stating  that  "the  two  oxime  groups  of  the  symmetrical¬ 
ly  structured  compound  are  identical  in  chemical  behavior  and  the  nucleo- 
pnilicity  of  their  anions  correspond  to  that  of  a  monoxime  with  pKa  =  7.80. 

The  fraction  of  the  oxime  anion  in  the  physiological  pH  range  is: 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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TOXOGONIN 


%  OXIME  ANION 

pH 

O 

O 

CXI 

37 "  ! 

M0N0XIME 

BIS0X IME 

SUM 

M0N0XIME 

BIS0X IME 

SUM 

6.80 

15.1 

3.1 

18.2 

24.0 

4.4 

28.4  | 

7.10 

26.2 

5.7 

31.9 

38.  7 

8.4 

47.  1  j 

7.20 

30.9 

7.1 

38.0 

44.3 

10.3 

54.6  j 

7.30 

36.0 

8.7 

44.7 

50.0 

12.6 

62.6 

7.35 

38.7 

9.7 

48.4 

52.9 

14.0 

66.9 

7.40 

41.5 

10.7 

52.2 

55.7 

15.4 

71.1 

7.45 

44.3 

11.9 

56.2 

57.4 

17.0 

74.4 

7.50 

47.1 

13.1 

60.2 

61.3 

18.6 

79.9 

7.60 

52.9 

16.0 

68.9 

66.6 

22.4 

89.0 

■SB 

58.5 

19.3 

77.8 

71  .5 

26.6 

98.1 

7.80 

64.0 

23.2 

87.2 

76.0 

31  .4 

107.4 

nr,  31 

There  had  been  a  report'’  that  TOXOGONIN  in  high  doses  was  capable  of 
overcoming  the  blood/brain  barrier  threshold  applicable  to  other  quaternary 
oximes. 

In  view  of  the  above  report3^  that  in  high  doses  TOXOGONIN  is  capable 
of  overcoming  the  blood/brain  threshold,  we  plan  some  further  lipophil icity 
experimentation  on  TOXOGONIN. 

However,  our  present  lipophil icity  experiments  reported  here  were 
carried  out  under  our  standard  concentration  conditions,  8x10-3  molar,  and 
10%  oil/water  solution.  These  drug  concentrations  and  a  10%  oil/water 
solution  have  proven  a  good  model  for  blood/brain  partitioning  of  drugs. 

These  standard  conditions  have  proven  reliable  to  allow  us  to  measure  the 
lipophil icities  of  more  than  50  other  diverse  compounds.  Our  standard 
conditions  have  allowed  us  to  measure  1 i poph ilicities  as  low  as  that  of  morphine 
over  the  range  from  pH  2  to  pH  1 1 ,  including  the  physiologic  range. 

Drug  Distribution  Coefficient 

Morphine  Ternp'  °7.10  °7.40  °7.70 

Sulfate  37°C  "OS'  HI  "OS" 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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As  a  reliable  model  for  penetration  of  a  drug  through  the  lipid  blood/ 
brain  barrier,  the  reliability  of  our  measured  morphine  1 i pophi 1 i ci ti es  was 
verified  by  admi ni steri ng  morphine  i.v.  totals  where  the  physiologic  pH's 
of  the  rats  were  varied  from  7.14  to  7.70.  Both  the  relative  ED50 1 s 

at  the  various  physiologic  pH's  and  the  relative  concentrations  of  radioactive 
morphine  in  various  portions  of  the  brain  as  a  function  of  physiologic  pH 
were  correctly  predicted  by  our  experimental  1 i pophi 1 i ci ty  values. 


The  ionization  constants  of  T0X0G0HIN  were  also  determined  at  30^0. 

T0X0G0NIN  does  not  seem  to  hydrolyze.  Retitration  of  the  compound 
with  HC1  retraces  exactly  the  titration  curve,  so  it  seems  that  the  oxime 
anion  in  alkalinic  media,  up  to  a  pH  of  11.50  and  over  a  time  interval  of 
^20  minutes,  does  not  undergo  an  irreversible  reaction. 

Since  there  are  also  two  overlapping  pKa's  in  T0X0G0NIN,  the 
computer  program  was  used  for  the  cal cul ati on  of  the  two  p Ka ' s .  The  results  at  30° 
are  as  follows: 


* 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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As  we  finished  our  pK  measurements  for  the  following  compounds:  2-PAM, 
HI-6,  TMB-4  and  Toxogonin,  at  three  different  temperatures  (20’C,  30 ' C  and 
37°C)  ,  we  calculated  the  temperature  dependences  of  their  pK  's.  The  pK 
of  many  non-carboxyl  i  c  organi  c  acids  are  temperature  sensitive  ffor  example^ 
phenol  becomes  stronger  by  0.012  units  pKa  for  each  degree  rise  in  temperature) 

The  results  are  the  following: 


Compound 

r  i 

2- PAM 

0.0123 

. .  '  ■  J 

| 

HI-6 

0.0131 

TMB-4 

0.0125 

0.0154 

TOXOGONIN 

0.0143 

0.0103 

See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures 
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We  then  began  systematically  to  study  the  class  of  methyl imidazol ium 
oxime  derivatives,  the  BHH  compounds;  in  numerical  order 

5.  BHH-1 

Microelectrometr ic  technique 


The  ionization  constants  determinations  were  made  at  37’C.  The  compound 
has  only  one  titratable  group. 

The  compound  did  not  hydrolyse  subjected  to  our  standard  conditions 
during  titration  and  retitration.  The  results  are  the  following: 

No  partitioning  occurs  in  the  oil  phase. 


Temperature 

Conditions 

PKa 

37°C 

Aqueous 

i  7. 76:0.01 

37°C 

1 0°5  oil /water 

|  7 .76:0.006 

We  calculated  the  percent  oxime  anion  at  various  pH's. 


PH 

j  Oxime  anion 

6.80 

9.88 

7.10 

17.95 

7.20 

21.59 

7.30 

25.75 

7.35 

28.00 

7.40 

30.39 

7.45 

32.87 

7.50 

35.46 

7.60 

40.89 

7.70 

46.55 

7.80 

52.30 

★ 


See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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6.  BHH-3 

Microelectrometric  technique 

The  shape  of  the  titration  curve  indicates  that  the  compound  has 
only  one  titratable  group.  The  retitration  curve  is  identical  to  the  titra¬ 
tion  one,  indicating  that  in  the  conditions  used  (C  =  8xlO"^M,  temp.  = 

37°C  and  time  -  20  minutes,  at  pH  values  up  to  11.70),  the  compound  does 
not  undergo  hydrolysis  and  the  formation  of  the  oxime  anion  is  a  reversible 
process.  The  compound  does  not  partition  in  the  oil  phase  under  our  stan¬ 
dard  conditions.  The  results  are  the  following: 


Temperature 

Condi tions 

pKa  | 

37°C 

7. 71  iO .01 

37°C 

7. 70i0.01 

_ l 

pH  j 

'!  Oxime  anion 

6.80  ! 

10.95 

7.10  i 

19.71  i 

7.20 

23.61 

7.30 

23.00 

7.35 

30.39 

7.40 

32.37 

7.45 

35.46  : 

7.50 

38.14 

7.60 

43.70 

7.70 

49.42 

7.80 

55.16 

See  Figure  1  of  this  report 


for  sketches  of  the  molecular  structures. 


Temperature 

Conditions 

PKa 

37°C 

Preliminary  measurements 
shift  w a 

Aqueous 

in  the  10%  oil/water 
s  negligible. 

7.82±0.01 

solution  for  the  pK 

_ 

PH 

%  Oxime  anion 

6.80 

8.72 

7.10 

16.00 

7.20 

19.35 

7.30 

23.19 

7.35 

25.31 

7.40 

27.55 

7.45 

29.90 

7.50 

32.37 

7.60 

37.60 

7.70 

43.13 

See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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8.  BHH-5 

Microelectrometric  technique 

This  compound  has  one  titratable  group.  It  appears  that  this  compound, 
due  probably  to  the  long  aliphatic  side  chain,  does  have  a  tendency  to 
partition  into  the  oil  phase.  Our  preliminary  measurements  indicate  that 
the  pKa  shift  (between  aqueous  and  oil/water  titrations)  is  still  quite 
small,  in  the  range  of  experimental  error.  Thus,  this  compound  is  a  prime 
candidate  for  a  spectrophotometric  measurement  of  partitioning. 


Temperature 

Conditions 

?Ka 

37°C 

Aqueous 

i  -  —  -  - - 

7 .88r0 . 006 

The  percent  of  oxime  anion  is: 


pH 

%  Oxime  Anion  ! 

6.30 

7.68 

7.10 

14.23 

7.20 

17.28 

7.30 

20.33 

7.35 

22.79 

7.40 

24.37 

7.45 

27.09 

7.50 

29.42 

7.60 

34.41 

7.70 

39.78 

7.80 

45.41 

See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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This  is  an  interesting  compound  since  it  is  the  only  oxime  that  we  have 
studied  to  date  that  has  a  large  partition  coefficient.  This  was  determined 

and  in  an 
are  as  fol 1 ows: 


from  the  pK  shift  (a  pK 

)  by  titration 

in  aqueous 

solution 

octanol/watir  mixture. 

aThe  values  for 

three  determinations 

SAMPLE  NO. 

ApKa 

P 

log  P 

1 

0.544 

25.0 

1.398 

2 

0.555 

25.9 

1.413 

3 

0.527 

23.7 

1.374 

Average 

Value 

0.542  ±  0.014 

r — i 

1-H 

-H 

CO 

CM 

1.395 

It  should  be  noted  that  we  are  reporting  a  partition  coefficient  (P) 
rather  than  a  distribution  coefficient  (D).  This  stems  from  the  fact  that 
to  derive  D  from  P  we  must  have  some  information  on  the  identity  of  the 
species  distributed  in  the  octanol. 

A  brief  comment  about  the  high  value  of  P  seems  in  order.  This  high 
value  probably  arises  from  the  fact  that  we  are  partitioning  into  octanol 
and  the  long  tail  attached  to  the  oxygen  is  an  octyl  group.  The  close 
structural  similarity  between  the  model  lipid  (octanol)  and  the  octyl  tail 
probably  accounts  for  the  high  P.  In  the  next  quarter,  experiments  using 
other  solvents  should  help  establish  whether  this  is  the  reason  for  the  high 
P  or  whether  other  factors  are  involved. 


*  See  Figure  1  of  this  report  for  sketches  of  the  molecular  structures. 
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9.  BHH-7 

Microelectrometric  technique 

This  compound  is  a  bis  pyridinium  derivative  (not  an  imidazi 1 ium  one) 
with  two  titratable  groups.  Color  changes  during  titration  frcm  colorless 
to  deep  yellow.  It  is  reversible  upon  retitration. 

The  pK  's  are  the  following: 

a 


Temperature 

Conditions 

pKa  j 

30°C 

30°C 

Aqueous 

Oi 1 /water 

GraDhical  Method  !  Computer  Method 

pKal 

pKa2 

pKal 

pKa2 

7.54-0.02 

7.48 

8.27i0.03 

8.28 

7 . 51 1 r0 . 01 

_ Z__ 

8 . 344-0 .01  ! 

pH 

7  Monoxime  Anion 

7  Bisoxime  Anion 

- 1 

Total  Oxime  Anion  j 

6.80 

16.23 

2.73 

19.06 

7.10 

27.96 

5.39 

33.35 

7.20 

32.32 

6.69 

39.52 

7.30 

38.09 

8.29 

46.33 

7.35 

40.34 

9.21 

50.05 

7.40 

43.65 

10.21 

53.36 

7.45 

46.49 

11.32 

57.81 

1  7.50 

49.37 

12.53 

61.90 

7.60 

55.11 

15.23 

70.39 

7.70 

60.71 

13.50 

79.21 

7.80 

66 .05 

22.22 

88.27  | 

E 


This  compound  has  two  titratable  groups:  one  attached  to  the  pyridinium 
ring,  the  other  to  the  imidazolium  ring.  During  titration  the  color  of  tne 
solution  turns  from  colorless  to  deep  yellow.  The  retitration  curve  retraces 
the  original  titration  curve,  so  no  hydrolysis  occurs  under  these  conditions. 


Temperature  Conditions  j 


Graphical  Method 


Computer  Method 


6.87  0.01  *7. 79  '  .005 : 6 .87 ' .C03i  7. 35'. 00 
6.35 


The  interesting  point  is  that  the  first  pKai  has  a  very  unusually  low 
value.  All  the  derivatives  studied  up  till  now  had  a  pKa-i  >  7.00  ana 

pKa2  >  8.10. 

It  would  be  interesting  to  establish  which  pKa  corresponds  to  which 
ring.  Single  oxime  imidazolium  ions  had  a  pKa  around  7.70  -  7.80.  One 
tends  to  assume  therefore  that  the  pKa  =  6.37  can  be  attributed  to  the  oxime 
group  attached  to  the  pyridinium  ring.  Although  unsubsti tuted  2-PAM  has 
a  pKa  of  7.68  versus  the  first  pi<a  3hH-3  has,  which  is  such  a  low  value, 
the  concentration  of  oxime  anions  of  3HH-3  is  very  high  at  physiological 
pH  levels: 


6.80 

7.10 

7.20 

7.30 

46 

63 

63 

73 

7.3 

7.4 

5 

0 

75 

77 

7.4 

5 

79 

7.5 

0 

81 

7.6 

0 

84 

7.7 

o 

87 

7.8 

0 

89 

54.34 
78.25 
36 . 56 


11.  BHH-1 3 


a.  Mi croel ectrometri c  technique 


This  compound  has  also  one  titratable  group--no  hydrolysis  occurs 
during  titration. 


pKal 


Conditions 


Aqueous 


7. 83*. 003 


T 
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b.  UV  spectrophotometr ic  technique 
The  character istic  optical  constants  tor  this  compound  ire: 


SOLUTION 

pH 

%  OXIME  ANION 

\  ( nm) 

max 1  ' 

r 

0 . IN  HC1 

1 

0 

111 

1.434  x  10  " 

0 . IN  NaOH 

13 

100 

292 

1.851  x  10H 

Phosphate 

Buffer 

7.4 

27.1 

278 

1.414  x  10  " 

The  distribution  coefficient  0  =  0.14  ±0.005.  This  distribution 
coefficient  is  based  on  the  concentration  of  the  oxime  in  both  the  aqueous 
and  in  octanol  phases.  If  we  had  only  used  the  measured  aqueous  concentra¬ 
tion  of  oxime  and  used  the  difference  between  this  value  and  the  total 
amount  of  oxime  used  as  the  octanol  concentration,  the  distribution  coeffi¬ 
cient  would  have  been  reported  as  0.24  ±  0.02. 


.a  (.k  t. t  t.a  l.l  ».l  l 


12.  BHH-14 


:NOH  •  W 


N  CH. 


ch2och2c-ch3 


a.  Microelectrometr ic  technique 

The  ionization  constant  of  this  compound  as  determined  by  electrometric 
titration  is 

pK  =  7.81  t  0.01  at  37°  C . 

a 


b.  UV  spectrophotometr ic  technique 

The  percentage  of  oxime  anions  in  the  physiological  pH  range  is  as  follows 

pH _ %  OXIME  ANION 

6.80  8.85 

7.10  16.22 

7.15  17.85 

7.20  19.56 

7.25  21.48 

7.30  23.48 

7.35  25.61 

7.40  28.01 

7.45  30.24 

7.50  32.72 

7.60  37.98 

7.70  43.53 

7.80  49.25 


•  *  f  •*'  f  •  _  a  '  —  •  4  a  (  .  >  _  *  ^  '  a  .  "a  ^  a  *”  a  **" „  •*” .  a**  a  *"  ■’»<*,  •*  ,  a’ 
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The  distribution  coefficient  of  this  compound  was  determined  by  UV 
spectrophotometry  since  the  partitioning  was  too  small  to  determine" 
electrometrically.  The  characteristic  optical  constants  are: 


SOLUTION 

pH 

l  OXIME  ANION 

A  ( nm) 

maxv  ' 

1 

0 . IN  HC1 

1 

0 

268 

1.190  x 

10" 

0 . IN  NaOH 

13 

100 

292 

1.473  x 

10" 

Phosphate 

Buffer 

7.4 

28 

274 

1.130  x 

10" 

The  distribution  coefficient,  D  =  0.012  ±  0.008. 


*  See  Figure  1  of  this  report  for  sketches  of  molecular  structures. 
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13.  BHH-15 


•  Cc" 


CH2OCHC-CH3 

ch3ch3 


a.  Microelectrometr ic  technic 


The  pK  value  as  determined  by  electrometr v  t-tration  is  7.78  t  0.02 

at  37° C.  a 

The  percentage  of  oxime  anion  in  the  physio’ogua’  range  is  as  follows 


PH 

%  OXIME  ANION 

6.80 

9.43 

7.10 

1  7 . 28 

7.15 

18.99 

7.20 

20.83 

7.25 

22.79 

7.30 

24.88 

7.35 

27.09 

7.40 

29.42 

7.45 

31.87 

7.50 

34.42 

7.60 

39.78 

7.70 

45.41 

7.80 

50.50 

b.  UV  spectrophotometr ic  technique 

The  distribution  coefficient  0  =  0.18  ±  0.02,  determined 
spectrophotometr ical ly  only.  The  optical  constants  are: 


SOLUTION 

pH 

%  OXIME  ANIONS 

X  ( nm) 

maxv  ’ 

Z 

0 . IN  HC 1 

1 

0 

269 

1.18  x 

10“ 

0 . IN  NaOH 

13 

100 

293 

1.47  x 

10" 

Phosphate 

Buffer 

7.4 

29 

275 

1.13  x 

10" 

*  See  Figure  1  of  this  report  for  sketches  of  molecular  structures. 
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VII.  Conclusions 

While  it  is  too  early  to  have  firm  conclusions  after  only  the  first 
year  of  a  several-year  new  research  project,  we  have  already  made  some 
significant  findings  which  lead  to  insight  into  the  characterization  of  the 
oximes  as  they  may  relate  to  their  in-vivo  efficacy  as  reactivators  of  AChE 
inhibited  by  organophosphorus  compounds. 

A.  General  conclusions 

Our  theoretical  and  experimental  results  to  date  indicate  that  we  are 
pursuing  fruitful  lines  of  investigation.  Our  results  to  date  already  shed 
considerable  insight  into: 

1.  Theoretical  stereoel ectronic  reequisites  for  an  effective  oxime 
reactivator  of  inhibited  AChE  (using  the  EMPC  map  of 
deprotonated  2-PAM  as  a  model); 

2.  The  fundamental  underlying  theoretical  basis  for  the  lack  of 

1 ipophi 1 ic ity  of  the  deprotonated  oximes  (even  though  they  are 
formally  overall  neutral  species),  since  the  EMPC  maps  indicate 
these  species  to  be  strongly  dipolar  ions; 

3.  Theoretical  techniques  to  predict,  prior  to  synthesis  of  a 
hypothetical  new  oxime,  whether  or  not  it  will  have  the  correct 
stereoelectronic  properties  to  reactivate  AChE  and  whether  the 
deprotonated  species  will  have  much  tendency  to  be  lipophilic 
and  partition  into  the  CNS  to  be  an  effective  reactivator  of 
inhibited  AChE  in  the  CNS; 

4.  A  strategy  using  the  principle  of  complementarity  of 
deprotonated  oximes  to  map  out  the  possible  stereoelectronic 
characteristics  of  the  AChE  active  site  (esteratic  and 
anionic).  This  is  significant  because  it  has  not  yet  been 
possible  to  crystallize  AChE  to  determine  experimentally  the 
structure  and  characterization  of  the  AChE  sites; 

5.  The  physicochemical  reason  for  the  non-ef fecti vi ty  of  the  oxime 
reactivators  in  the  CNS  (since  neither  the  quaternary  nor  the 
deprotonated  species  are  very  lipophilic  and  thus  will  not  have 
much  tendency  to  partition  through  the  CNS  blood-brain  barrier 
under  most  physiological  conditions). 

B.  Specific  conclusions 

Deprotonated  2-PAM  (the  oximate  that  has  been  suggested  to  be  the 
pertinent  species)  is  formally  an  overall  neutral  species  [with  a  formal  + 
charge  from  the  quaternization  and  a  formal  -  charge  from  the  deprotonati on 
of  the  oximej . 

However,  the  three  dimensional  isopotential  EMPC  map  around 
deprotonated  ( zwi tter ion  ic )  2-PAM  (generated  from  our  calculated  quantum 
chemical  wave  function  of  deprotonated  2-PAM)  indicated  that  dearly  there 
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is  a  strongly  positive  region  and  a  strongly  negative  region.  Thus,  this 
species  appears  to  be  a  strongly  dipolar  ion. 

This  has  significant  implications  for  several  aspects  of  the  use  of 
oxime  reactivators.  First,  the  size  and  shape  and  relative  spatial 
relationship  of  the  positive  and  negative  isopotential  portions  of  these 
EMPC  maps  should  have  a  complementary  relation  to  the  spatial  relationship 
and  electronic  character istics  of  esteratic  and  anionic  subsites  of  the 
active  site  of  AChE.  The  is  particularly  important  in  view  of  the  complete 
lack  of  molecular  crystal  structural  information  on  any  AChE.  Our  isopoten¬ 
tial  EMPC  maps  of  deprotonated  2-PAM  can  help  shed  insight  onto  the  size, 
shape  and  stereoelectronic  characteristics  of  the  AChE  active  sites-- 
quantities  at  present  not  yet  obtained,  since  so  far  it  has  apparently  not 
been  possible  to  crystallize  AChE  to  determine  its  crystal  structure. 

Second,  the  isopotential  EMPC  map  of  deprotonated  2-PAM  is  as  unique  as 
a  fingerprint  in  identifying  stereoelectronic  character  istics  that  have  the 
proper  stereoelectronic  complementarity  and  are  favorable  for  effective 
reactivation  of  AChE  inhibited  by  organophosphorus  compounds.  The  cor¬ 
responding  calculated  isopotential  EMPC  maps  of  the  oximates  of  proposed 
hypothetical  new  oxime  reactivators  can  be  computer- compared  to  those  of 
deprotonated  2-PAM  prior  to  the  synthesis  of  the  new  compounds  to  see  if 
these  new  compounds  have  the  proper  stereoelectronic  characteristics  to  be 
effective  new  reactivators . 

Third,  our  three-dimensional  isopotential  EMPC  maps  around  deprotonated 
2-PAM  also  shed  insight  into  another  important  aspect  of  the  properties  of 
2-PAM  and  deprotonated  2-PAM.  2-PAM  itself  is  a  positively  charged  quater¬ 
nary  species.  As  such,  it  is  not  expected  to  be,  nor  is  it,  very 
lipophilic.  2-PAM  has  little  tendency  to  partition  from  aqueous  to  lipid. 

It  had  been  hoped  that  when  2-PAM  was  deprotonated  and  was  then  an 
overall  neutral  species,  that  it  would  have  a  much  greater  tendency  to 
partition  into  lipid. 

Our  experimental  determination  of  the  1 ipophi 1 icity  of  2-PAM  as  a 
function  of  pH  indicated  that  even  when  the  positively  charged  quaternary  2- 
PAM  was  deprotonated ,  the  resulting  species  was  still  only  very  slightly 
lipophilic.  This  behavior  is  commensurate  with  our  EMPC  map  results, 
namely,  that  deprotonated  2-PAM  (while  overall  neutral)  has  the  electronic 
characteristics  of  a  strongly  dipolar  ion. 

Because  2-PAM  would  have  little  tendency  to  partition  through  the 
lipophilic  blood-brain  barrier  under  usual  physiological  conditions,  2-PAM 
is  not  an  effective  reactivator  in  the  CNS  of  AChE  inhibited  by  organophos¬ 
phorus  compounds. 

Our  EMPC  map  of  deprotonated  2-PAM  indicated  that  while  this  species  is 
neutral  overall,  it  appears  to  be  a  strongly  dipolar  ion  and  thus  would  not 
be  expected  to  have  a  significant  tendency  to  partition  into  lipids. 

Thus,  our  theoretical  finding  that  the  EMPC  map  of  deprotonated  2-PAM 
indicates  it  to  be  a  strongly  dipolar  ion  has  concomitantly  been  validated 
by  our  experimental  1 ipophi 1 icity  determination  of  2-PAM  as  a  function  of 


pH.  Neither  the  quaternary  2-PAM  nor  the  deprotonated  2-PAM  has  much  ten¬ 
dency  to  partition  from  aqueous  to  lipid.  Thus,  under  most  physiological 
circumstances,  2-PAM  would  not  be  expected  to  partition  into  the  CNS. 

One  of  the  Army's  interests  is  in  finding  and  developing  new  oxime 
reactivators  that  would  be  effective  in  the  CNS.  Oximes  that  were 
lipophilic  would  have  a  better  chance  to  be  effective  reactivators  in  the 

CNS. 

If  the  EMPC  map  of  the  deprotonated  species  of  a  proposed  new  oxime 
reactivator  indicated  that,  like  deprotonated  2-PAM,  it  was  a  strongly 
dipolar  ion,  then  this  oxime  would  not  appear  to  be  a  good  candidate  for  an 
effective  reactivator  of  inhibited  AChE  in  the  CNS.  (While  there  have  been 
reports  that  the  CNS  blood-brain  barrier  loses  some  of  its  integrity  after 
exposure  to  organophosphorus  AChE  inhibitors,  sufficient  blood-brain  barrier 
integrity  must  still  remain,  because  2-PAM  has  been  shown  not  to  be  an 
effective  AChE  reactivator  in  the  CNS.) 

We  are  continuing  our  theoretical  and  experimental  investigations  of 
other  oxime  reacti vators . 

Our  experimental  microelectrometr ic  results  on  the  pH  dependence  of  the 
1 ipophi 1 icities  of  the  quaternary  oximes  and  their  deprotonated  species 
indicate  that  neither  the  positively  charged  quaternary  nor  the  overall 
neutral  deprotonated  species  are  very  lipophilic.  However,  in  view  of  the 
fact  that  inhibited  AChE  can  be  reactivated  by  low  concentrations  of  the 
oximes,  we  set  up  and  are  using  a  more  sensitive  UV  spectrophotometric 
method  in  which  the  two  phases  are  separated  by  centrifugaaaon  and  the 
amount  of  material  in  each  phase  is  determined  from  calibration  curves. 
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VIII.  Recommendations 


General 


To  continue  to  mesh  together  our  theoretical  quantum  chemical  and  EMPC 
map  results  with  our  experimental  1 ipophi 1 ic i ty  results  to  get  a  better 
understanding  of  the  efficacy  of  the  oxime  reacti vators . 

B.  Theoretical 

Our  recommendations  are: 

1.  Pyridinium  oximes  and  bis-pyridinium  oximes 

To  carry  out  in  depth  further  quantum  chemical  calculations  for  the 
conformational  energy  profile  of  2-PAM  as  a  function  of  various  geometry 
variations  based  both  on  the  new  experimental  crystal  structure  of  2-PAM  and 
on  our  ansatzed  structure,  and  to  compare  the  results  carefully.  Since 
there  are  still  very  few  experimental  crystal  structures  of  these  quaternary 
oxime  reactivators,  this  is  an  important  point.  There  are  no  experimental 
crystal  structures  for  most  of  the  oximes  of  interest,  including  the  BHH 
compounds.  Thus,  for  quantum  chemical  calculations  on  the  newer  oximes,  the 
structures  will  have  to  be  ansatzed  as  a  starting  point  for  further  quantum 
chemical  calculations. 

To  carry  out  further  in  depth  quantum  chemical  calculations  on 
deprotonated  2-PAM  (the  oximate)  using  the  new  crystal  structure  of  2-PAM  as 
a  starting  point  and  making  extensive  geometry  variations,  including  bond 
lengths  and  bond  angles,  to  calculate  the  conformational  energy  profile. 

The  covergence  of  the  quantum  chemical  wave  function  of  this  species  has 
proven  to  be  sensitively  dependent  on  the  geometry.  The  oximate  of  2-PAM  is 
the  model  for  all  other  effective  oxime  react i vators .  Thus,  a  thorough 
understanding  of  its  wave  functions  and  behavior  will  be  an  important  guide 
for  the  oximates  of  the  other  oximes. 

For  both  2-PAM  and  its  oximates,  the  EMPC  maps  will  be  generated  for 
various  conformations  to  assess  their  sensitivity  to  geometry. 

To  carry  out  in  depth  quantum  chemical  investigations  of  the  conforma¬ 
tional  energy  profile  of  HI-6  as  a  function  of  various  geometry  variations, 
to  calculate  conformational  energy  profiles,  based  both  on  the  new  ex¬ 
perimental  crystal  structure  and  on  our  ansatzed  geometry,  and  to  compare 
the  results  carefully.  HI-6  is  the  prototype  of  the  bis-pyridinium  oximes 
which  have  many  more  conformational  degrees  of  freedom  than  the 
monopyridini urn  compounds. 

To  carry  out  similar  quantum  chemical  calculations  for  the  conforma¬ 
tional  energy  profiles  of  deprotonated  HI-6  (the  oximate)  as  a  function  of 
geometry  variations  based  both  on  the  new  crystal  structure  and  the  ansatzed 
structure,  and  to  compare  the  results  carefully.  These  are  models  for  the 
bis-pyridinium  compounds  which  have  many  more  degrees  of  conformational 
freedom  than  the  monopyridinium  compounds. 


Then,  to  generate  EMPC  maps  around  various  conformations  of  HI-6  and 
its  oximates  and  to  study  the  variations  of  the  EMPC  maps  as  a  function  of 
conformation. 

2.  Interactions  with  AChE 

To  use  the  principle  of  complementarity  of  the  EMPC  maps  around 
deprotonated  2-PAM  (the  oximate)  to  infer  stereoel ectronic  characteristics 
of  the  AChE  active  site  (esteratic  and  anionic)  which  have  not  yet  proven 
amenable  to  experiment. 

To  carry  out  quantum  chemical  calculations  and  intermolecular  interac¬ 
tion  calculations  of  deprotonated  2-PAM  interacting  with  a  model  of  the  AChE 
active  site. 

To  carry  out  similar  calculations  of  deprotonated  2-PAM  interacting 
with  inhibitng  AChE. 

3.  To  formulate  how  the  quantum  chemical  and  EMPC  map  information 
would  be  used  to  set  up  computer  based  structure-activity 
prediction  program  for  the  efficacy  of  oximes. 

This  will  be  done  based  on  EMPC  maps  as  we  get  theoretical  results  on 
more  molecules. 

4.  To  use  the  quantum  chemical  and  EMPC  map  information  to  set  up 
a  program  for  lead-directed  synthesis  of  new,  more  effective 
oximes. 

We  will  also  be  doing  this  based  on  the  EMPC  maps  of  already  effective 
oximes  and  how  they  vary  with  reactivating  effecti vities.  Since  quantum 
chemical  calculations  can  be  run  and  EMPC  maps  generated  for  proposed 
hypothetical  oximes  prior  to  their  synthesis,  this  theoretical  approach  is  a 
very  effective  predictive  tool. 

5.  BHH  Compounds  -  quaternary  methyl  imidazolium  oximes 

To  carry  out  the  quantum  chemical  calculations  and  EMPC  map  generation 
on  these  compounds  as  time  and  resources  permit.  In  particular,  to  compare 
these  to  2-PAM  and  its  oximate  as  a  model  effective  oxime  reactivator  and  to 
compare  the  BHH  effective  reactivators  against  those  which  are  "protective". 

C.  Experimental 

To  continue  to  carry  out  the  determinations  of  1 ipophi 1 1 i ci ty  of  oximes 
by  both  the  microelectrometr ic  and  the  UV  spectrophotometric  techniques  as  a 
function  of  pH. 

Further,  for  the  one  species  to  date,  BHH-5,  that  had  a  reasonable 
1 ipophi 1 icity  in  octanol,  to  attempt  to  identify  the  specific  species  which 
have  partitioned  into  octanol  using  NMR  technique. 
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